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I. INTRODUCTION 

The problem of quantum state transfer along the spin chain [l| acquires significance in the 
study of quantum communication systems. Several models of the end-to-end quantum state 
transfer along the spin-1/2 chain were suggested |2|-l9|. All these models are based on using 
either inhomogeneous chains with different coupling constants between nearest neighbors 
(NN) or homogeneous chains with different Larmor frequencies on the chain nodes. Seminal 



works l2- 



m^ 



61] use NN interactions in the inhomogeneous symmetrical chain. All coupling 
constants of such chains have fixed values, allowing one to transfer the state duringthe time 
interval independent of the length of the chain. The method proposed in refs.p, ISJ and 
developed in refs. [T], l8| is based on the week end bonds in the chain when interactions of all 
remote spins are taken into account (the week end bonds method (WEBM)). Transfer time 
increases with the length of the chain in this case. 

Since it is very difficult to construct inhomogeneous chains, an alternative method of 
state transfer along the spin chain has been proposed recently |9|. The effect of the high 
probability state transfer (HPST) along the homogeneous spin-1/2 chain is achieved here 
because of the specially tailored external inhomogeneous magnetic field. In the simplest 
case one needs only two nonzero end Larmor frequencies and zero inner ones (end Larmor 
frequencies method (ELFM)). 

The problem of quantum state transfer requires solving some physical questions. In 
particular, it is not clear whether the presence of inner nodes (so-called "body" of the chain) 
improves parameters of state transfer (transfer time and probability). How will the transfer 
time be changed if we fix end nodes and remove the body? Even this very simple question 
requires a special study and will be considered for two methods of state transfer (WEBM 
and ELFM) applied to the spin-1/2 chain governed by either an XY or XXZ Hamiltonian 
with dipole-dipole interactions. 

In this article, we develop an approach to the exact diagonalization of the XY-Hamiltonian 
with NN interactions for the combination of WEBM and ELFM, FiglH We show also the 
formal equivalence of the spin dynamics governed by the XY and XXZ Hamiltonians in the 
specially tailored inhomogeneous external magnetic field in the problem of quantum state 
transfer with a single excited spin. 

We perform the numerical simulations of quantum state transfer along the ten-node 
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FIG. 1: Two methods of the state transfer along the spin-1/2 chain; a is the distance between 
nearest neighbors in the body of the chain: (a) WEBM, ah Larmor frequencies are zeros; {h) 
ELFM along the homogeneous chain; the end node Larmor frequencies equal to O, while inner 
ones are zeros. 

chain with dipole-dipole interactions governed by different Hamiltonians (XY or XXZ) using 
different methods of state transfer (WEBM or ELFM) and different types of interactions 
(NN or all node interactions). We concentrate on the following problems in numerical 
simulations: 

1. Study the influence of the inner nodes on the transfer time over the given distance 
L. We will show that inner nodes may either descrease or increase the transfer time 
in dependence on the Hamiltonian, the method of the state transfer and the type of 
interactions. 

2. Compare the end-to-end transfer times along the chains governed by different Hamil- 
tonians with different methods of the state transfer and different types of interactions. 

Thus, the structure of this paper is following. We give some details regarding the spin 
dynamics in SecHIl The analytical description of the spin dynamics along the chain governed 
by the XY Hamiltonian with NN interactions using the combination of WEBM and ELFM 
is given in Sec. Illli The numerical simulations of the quantum state transfer correspond- 
ing to different Hamiltonians, methods of the state transfer and types of interactions are 
represented in SecJIVi We briefly summarize our results in the concluding section, SecJVl 



II. SINGLE EXCITED SPIN DYNAMICS 

Consider the spin dynamics along the spin-1/2 chains governed by either XY or XXZ 
Hamiltonian in the inhomogeneous external magnetic field. These Hamiltonians read: 

N ^ 2fc 

T-lxY = y , Dij{li,xlj,x + k,ylj,y) + / , ^i h,z, Dij = —^, (l) 

N N 

'HxXZ = 2_^ ^i,ji^i,x^j,x + h,ylj,y — '^h,zlj,z) + 2_^^i "^*'-^' (^) 

Here Dij is the coupling constant between the ith and the jth nodes. Hereafter we will 
study the one-dimensional spin-1/2 chains and use the dimensionless time r, distances ^„,m, 
coupling constants (i„,„ and Larmor frequencies a;„, defined as follows: 

r D 1 

_ n -A c — "'"* J — "i"* _ ^ _ 1 ^Q^ 

^^1,2 ^1,2 ^^, 



,,m 



, XY _ "n , ,XXZ _ "n 

-^^1,2 -^^1,2 

Using definitions (E]), the Hamiltonians (II]) and ([2]) may be written as follows: 

TV N 

TixY = DiHxYi y-XY = 2_^ dj^i{li^xlj,x + h,ylj,y) + /_^^i h,z, (4) 

».J=i i=l 

N N 

T^XXZ = D{HxXZi T^XXZ = ^^ di,j{h,xlj,x + h,ylj,y — '^Ii,zlj,z) + ^^i^j -^i,z- (5) 






In order to characterize the effectiveness of the quantum state transfer from the first to 
the iVth node, the fidelity F{t) has been introduced [l|: 

where fn,m is a transfer amplitude: 

N 
Jn,m\T) \Tn\C |?2y y^UnjUfnjC , Jn,m Jm,n- \' ) 

J=l 

Here Mj^-, i,j,= 1, . . . , iV, are components of the normalized eigenvector Uj corresponding 
to the eigenvalue \j of the matrix D, H = ^D: Duj = XjUj, where H is the matrix 
representation of Ti. 



A. Equivalence of XY and XXZ Hamiltonians 

It is simple to show the equivalence of the single excited spin dynamics governed by the 
XY and XXZ Hamiltonians. In fact, two Hamiltonians (jl]) and ([5]) differ by the diagonal 
parts, i.e. 

N N 

An = UxY - UxYz = 5Z(wf "^ - ^f''^)kz + Y. '^'^^^okzh,z- (8) 

Since we consider the single excited state transfer and both Hamiltonians commutes with 
Iz {z projection of the total spin), the spin dynamics is described by only one block of 
the Hamiltonian which may be written on the basis of A^ eigenvectors |ra), n = 1, . . . ,N, 
where notation \n) means that the nth node is excited, i.e. directe d op posite to the external 
magnetic field. Then the matrix representation AH of AT-t reads j 
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11| as follows: 



A^ = w^^ - w^^^ -A + TEn, (9) 

cj^^ = diag(a;f ^, . . . , w^^), w^^^ = diag(a;f ^^, . . . , cj^^^), A = diag(Ai, . . . , A^), 

TV ^ N N 

a^=j: ci,„, t=iy: dM + \ Y.^^r' - ^n- 

Thus both Hamiltonians are equal to each other if AH = 0. Since the number of the diag- 
onal elements in AH coincides with the number of Larmor frequencies in each Hamiltonian, 
this condition may be satisfied by the proper choice of the Larmor frequencies in either the 
XY or XXZ Hamiltonian. This conclusion remains correct for an arbitrary spin-1/2 system. 
Nevertheless, we consider XY and XXZ Hamiltonians independently as far as there is no a 
simple method of generation of the inhomogeneous magnetic field. 

Remember that F may be put to zero by the proper choice of the constant homogeneous 
magnetic field jj. For this reason, hereafter we consider the probability of the state transfer, 
^i'^) — l/iA^('r)P; instead of the fidelity as a characteristic of the state transfer effectiveness. 

III. THE EXPLICIT SOLUTION FOR THE SPIN DYNAMICS IN THE CHAIN 
GOVERNED BY THE XY-HAMILTONIAN WITH NN INTERACTIONS 



In this section, we develop the ideas of Ref. IQl] and give an analytical description of 



the single excited spin dynamics in the chain governed by the XY Hamiltonian with NN 



d, 



i,i+l 



interactions combining both WEBM and ELFM. In other words, one has to substitute 

6, l<i<N-l ^y f 0, 1 < i < A^ 

1, i = l, N-1 ^u, i = l, N 

into the Hamiltonian (jlj), which reads: 

TixY = 

N-2 ^ 

j=2 



(10) 



(11) 



Since we are interested in the single excited spin dynamics, the Hamiltonian has A^xA^ matrix 
representation (see SecHIl)), which reads (up to the scalar term, which is not important in 
our case and may be removed by adding a proper homogeneous constant external magnetic 
field): 

/2a; 1 ■■■ o\ 
1 S ■■■ 
1 5 ■•■ 



H 



D, D 



(12) 



■■■ 1 
\0 ■■■ 1 2a;/ 

First, one has to find eigenvalues and eigenvectors of D. For this purpose, we solve the 
equation 

Duj = XjUj, Uj = (uij, . . . , unj) (13) 

for the components of the eigenvectors Uj, where Xj are the solutions to the characteristic 
equation 

det{D - XEn) = 0, (14) 

and En is an N x N unit matrix. The system ( !T3|) may be written as follows: 

2uJUij + U2j = XjUij, (15) 

Uij + 5U3J = XjU2j, (16) 

6u(^k-i)j + Su(^k+i)j = XjUkj, 3 < k < N -2 (17) 

SU{N-2)j + UNj = XjU{N-l)j, (18) 

U{N-l)j + 2uUnj = XjUNj- (19) 



Eqs.f ll7p are known to have the following solution 



Ukj = Cxje-'^"' + C2Je'^P^ 2 < fc < AT - 1, 
Aj = 2(5 cos pj. 



(20) 
(21) 



Then eq. (fT^ becomes equivalent to the compatibility condition of the system of four equa- 
tions, eqs. (I15|16|18|IT9|1 for any j. After some transformations, this compatibility condition 
reads 



4z I o(2a; — 2o cosp) cos — -^ — p + cos — - — p ) x 

Ar-3 



(22) 



A^- 1 
5(20; — 25 cosp) sin — -^ — "p + sin 



p =0. 



2 ' 2 

One can show that this equation can be treated as the A^ degree polynomial equation for 
the variable X = cosp, and consequently, it has A^ roots. All roots can be separated into 
two sets of A^^i and N — Ni roots, which are solutions to one of two following equations: 



iV- 1 A^-3 

S(2uj — 26 cospj) cos — - — pj + cos — - — pj = 0, j = 1, . . . , A^i, Ni = N — 



, N-1 
o{2u — 26 cospj) sin — - — pj + sin 



2 

A^-3 



N 
~2 



(23) 



p, = 0, j = Ari + l,...,Ar. 



2 '' 2 

Here [a] means the integer part of a. Thus, using eqs.( l23l) . expressions for the normalized 
eigenvector components can be written as follows: 



M,co.<^-^'P'., = l 



u 



Ij 



6A^ sin 



Aj cos 



Uki 



UNj 



Aj sin 



(N - l)p, 

2 
^^ {N + l-2k)pj 



j=N, + l,...,N 



(24) 



(A^ + 1 - 2k)pi , 



M,cos(^^, , = 1 
. {N-l)p, 



, k = 2,...N~l, j = N, + l, 

,iVi 



A^ 



-6A^ sin ^^——^, j = Ni + 1,...,N 



where 



A = < 



+ 5^(1 + cos(Ar - l)pj) + \ ^ ^^^ ) , J = 1, 



2 

A^-2 



6\1 - cos{N - l)pj) - 



2sinj9j 

sin(A^ - 2)pj 

2 sin Pj 



-1/2 



iVi 



.(25) 



j = Ari + l,...,Ar 



Then the probabihty of the end-to-end excited state transfer reads: 

2 



P{t) = (Ar|e-*^^^^|1) 



N 






-iTXj/2 



(26) 



Ni 



J^A^jCOS^N - 1)^6-'^^''°'' 



N 



\ rj ^—iSr cos pj 



j=l " j=Ni+l 

which will be used in those examples of Sec JIVI where NN interactions are considered. 



A. Spin chain of four nodes 

The derived formulas become most simple in the case A^ = 4. The eigenvalues (I2T]) have 
the following explicit forms: 



_ 2UJ-6+ ^{2uj + 6y + 4 _2uj-6~ ^(2tj + (5)2 + 4 



(27) 



_ 2UJ + 6+ v/(2w-5)2 + 4 _2uj + 6- y/{2uj-6y + 4: 
^3 — ;:; 5 -^4 — • 



Formula (126|) reduces to the following one: 

1 



Pir) 



-e-^'^l^ (cos(r/3(a;, 5)/4) - zf^^ sin(r/3(a;, 5)/4) 
2 [ p(w,<^) 

-le^-^-/^ |cos(™(a;, 5)/4) - i^^^r^ sin(ra(a;, 5)/4)^ 



(28) 



where 



a 



a{uj, S) 



(w, S) = v/(2a; + 5)2 + 4, /3(a;, 5) = ^(20; - 5)2 + 4. 



(29) 



Let us demonstrate that the perfect state transfer is possible at some time moment tq if the 
parameters u and 6 have been properly chosen. First, we remark that perfect state transfer 
is achieved if 



or 



cos(ro(5/4)| = l,cos(roa(w,(5)/4) = ±1, cos(ro/3(a;, (5)/4) = =f1, 



sin(ro(5/4)| = 1, cos(roa(w, 5)/4) = ±1, cos(ro/3(a;, (5)/4) = ±1. 



(30) 



(31) 



Eqs.( l30|) and ( 1311) are equivalent to the following pair of complete systems of algebraic 
equations (remember that tq, 6, a, and (3 must be positive, while u must be real): 



tqS — Anrii, Toa{co, 6) = 47rn2, rof3{u, 6) = Ann^ 



(32) 



ni = 1, 2..., n2 = 1, 2, ..., n2 - n^ = ±1, ±3, ±5, ..., 
or 

tq5 = 7r(4ni + 2), Toa{u, 5) = 47rn2, to/3{u, 5) = Attus, (33) 

ni = 0, 1, ..., ns = 1, 2..., na - ns = 0, ±2, ±4.... 

If w = (WEBM), then we obtain alternating spin chain of four nodes. In this case, 
a{0,6) = f3{0,6) = v^M-4 so that system ( l32l) becomes inconsistent, while system ( 133|1 
acquires the following form: 

To6 = 7r(4ni + 2), roa(0, 6) = toVS^ + 4 = 47rn2, (34) 

ni = 0,1,..., n2-ni = 1,2..., 



which agrees with the results obtained for the alternating spin-1/2 chain of four nodes 10 |. 
The relation between rii and n2 in eq.l lMj) is a consequence of the inequality 6 < a(0, 6). 

IV. NUMERICAL SIMULATIONS OF THE SPIN DYNAMICS. 

Now we represent results of the numerical simulations of the end-to-end single excited 
state transfer along the spin-1/2 chain of ten nodes. Considering the problem of state 
transfer between two end nodes separated by the dimensionless distance L = C,i,n, first, we 
recall the simplest model of two spins with the distance L between them (Sec JIVv^jl . Next, 



in Sec JIVBl we find end-to-end transfer times over the same distance L along the ten-node 
chain using both XY and XXZ Hamiltonians, with both WEBM and ELFM using both NN 
and all node interactions. In particular, we study the problem whether the presence of inner 
nodes reduces the time interval required for state transfer over the distance L in comparison 
with the above simplest system of two interacting spins. It will be shown that the answer is 
not always positive. 

Note that we solve the optimization problem for the case of all node interactions (i.e., we 
find such parameters of the spin chain which provide the end-to-end HPST during as short 
as possible a time interval). Then, we consider the spin dynamics using NN interactions and 
keeping the same parameters of the spin chain. This allows us to see whether NN interactions 
give different transfer times in comparison with the case of all node interactions. 



A. The state transfer over the distance L 

n 

It is well known [2] that the quantum state can be perfectly transferred between two 
spin-1/2 nodes separated by the distance L = ^i^v during the dimensionless time interval T2 
such that (see eq.Q) 

T,{L) = -^ = L\. (35) 

However the time interval T2 increases with an increase in L so that the direct node-to-node 
quantum state transfer over the long distance L becomes impossible because of the quantum 
decoherence. It is assumed that the set of inner nodes placed between the above two nodes 
may help one to overcome the problem of decoherence. However, it seemed that these nodes 
were not always helpful. 

We will use the dimensionless time interval T2{L) as a characteristics of the A^-node chain. 
We say that the inner nodes in the A^-node chain are usefull if the end-to-end transfer time 
along this chain is less than T2. 

It is remarkable that the parameter T2{L) depends only on L and does not depend on the 
Hamiltonians and methods of state transfer. In fact, the XY Hamiltonian coincides with 
the XXZ Hamiltonian up to the scalar term in the matrix representation (see eq.Q with 
Qf^ = Q^"^ , Qf^^ = ^^^^ , N = 2) a.s far as the problem of the single quantum state 
transfer along the two-node chain is considered. This term does not effect on the probability 
of state transfer. 

B. The ten- node chain 

We represent results of the numerical simulations of the end-to-end state transfer along the 
ten-node chain using both WEBM and ELFM. Considering WEBM, we will use parameter 
6 introduced in eq.( 1T0|) and zero Larmor frequencies. Choosing the value of 6, we follow the 
ideology of ref . [7| , where the distance between the end nodes and the body of the chain is 
twice as long as the distance between neighbors in the body. Thus we take 6 = 8, see FigJT^ 
with a = 1. However, we have observed that all conclusions remain valid for different values 
of 6 as well. Considering ELFM, we take 6 = 1 and nonzero end Larmor frequencies u, while 
all inner Larmor frequencies are zeros, see FigHb with a = 1 and Q = u. 
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Results of our numerical simulations are represented in Figsl2}l5l namely, the functions 
P{t) obtained using WEBM for the spin-1/2 chains governed by the XY and XXZ Hamil- 
tonians are shown in FigsjJ] and [3] respectively; similarly, functions P{t) obtained using the 
ELFM for the spin chains governed by the XY and XXZ Hamiltonians are shown in FigsH] 
and [5l respectively. Dynamics with both all node interactions and NN interactions are rep- 
resented therein. Parameters T^'^"), P^'*") = P(T('^")) and T^™), p("") = P{T^'''')) are the 
end-to-end transfer times and probabilities of these transfers in the case of all node interac- 
tions and NN interactions respectively. In the following, we collect the results of analysis of 
all figures EB 

(i) Let T^WEBM) ^^^ rpiELFM) ^^ ^Yie transfer times corresponding to WEBM and ELFM 
respectively. Comparing the transfer times shown in Figsl2}l5]with appropriate interval T2 we 
conclude that the inner nodes do not always reduce the transfer time in comparison with T2 
(see Table I). We see that WEBM with an XY Hamiltonian is the most promising method 





all node interactions 


NN interactions 


XY 


XXZ 


XY 


XXZ 


rp{WEBM) 

r2(ll/2) 

rp(ELFM) 

T2(9) 


0.041 
0.319 


1.262 

0.144 


0.055 
212.017 


272.228 
21.194 



TABLE I: Effect of the inner nodes on the transfer time 



of state transfer, namely, such chains allow one to transfer the state over the fixed distance 
L (in this case, ^1,2 = 6,10 = 1, 6,i+i = 1/2, 2 < i < 8, so that L = X]Li^i,«+i = 11/2) 
during the shortest time interval ^"^^""^^ = 0.041 see Fig|2^ and Table I. However, one has to 
remember that the XXZ Hamiltonian is the most natural one describing the spin dynamics 



in the strong external magnetic field. ELFM is better suited for this Hamiltonian. It gives 
^-^ = 0.144 (in this case, ^j,j+i = 1, 1 < i < 9, so that L = J2i=i 6,j+i = 9), see Figj5|i and 
Table I. 

(ii) Comparing T*^""^ with T*^"") (which are shown in Figsj2|L-[5tL andlJb-IHb, respectively) we 
see that T*^""^ > T*^"") in all considered experiments with dipole-dipole interactions, which 
is reflected in Table II. This observation may be explained as follows. As we have seen, the 
probability of the end-to-end state transfer P{t) = |/iAr(r)p (where /itv is defined by eq. ([7])) 
is a superposition of the oscillating functions. Figsj2}l5] demonstrate us that the oscillation 
with the minimal frequency has the maximal amplitude. The transfer time is defined mainly 
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WEBM 


ELFM 


XY 


1.334 


664.444 


XXZ 


215.710 


146.929 



TABLE II: Table of values 



n~'{an) 



by this oscillation. Note that this frequency is simply related with the minimal eigenvalue 
Amin of the matrix D (see eq.( !T2l) ) in the case of XY Hamiltonian with WEBM and NN 

= 26.648, while 



interactions UM- In fact, Xmin = 0.118 in this case, so that T*^"") ^ y^^~ 

the calculated value of T*^"") is 28.698, see Figj2b. Remote node interactions increase this 

minimal frequency, which leads to the decrease in transfer time. 

(iii) Let T^^"^^ and T^^^^'^ be the transfer times corresponding to the chains governed by 

the XY and XXZ Hamiltonians respectively. Comparing the transfer times in Fig|2] with 

the appropriate times in FigJ3]and the transfer times in FigJUwith the appropriate times in 

Figjniwe conclude that 



WEBM : T^^^) < T^^^^^, 
ELFM: T^xy)>T'^xxz) 



(36) 



see Table II L 





all node interactions 


NN interactions 


WEBM 

ELFM 


0.033 

2.212 


2.0 X 10^^ 
10.004 



TABLE III: Table of values 



rp{XXZ) 



(iv) Comparing the transfer times in Figs J2][5] with the appropriate times in FigsJUS] we 
conclude that WEBM is preferable for the chains governed by the XY Hamiltonian, while 
ELFM is suitable for the XXZ Hamiltonian, see Table IV. Since we use different L for 
WEBM (L = 11/2) and for ELFM (L = 9), we have to introduce parameter b = (18/11)3 
in order to get the correct comparison. 
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FIG. 2: The end-to-end HPST along the chain often nodes governed by the XY Haniiltonian using 
WEBM; 6 = 8, r2(ll/2) ^ 522.682; (a) ah node interactions, T^"") = 21.518 < T2, P^"") = 0.976; 
(b) NN interactions, T(™) = 28.698 < T2, P^"") = 0.972 




FIG. 3: The end-to-end HPST along the chain often nodes governed by the XXZ Hamiltonian using 
WEBM; 6 = 8, r2(ll/2) ^ 522.682; (a) all node interactions, T^'*") = 659.630 > T2, P("") = 0.995; 
(b) NN interactions, T^™) = 142288.896 > T2, P^"") = 1.000 
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FIG. 4: The end-to-end HPST along the chain often nodes governed by the XY Hamiltonian using 
ELFM; uj = 2.203, T2(9) ^ 2290.221; (a) ah node interactions, T^"") = 730.786 < rs, P^"") = 0.985; 
(b) NN interactions, T^™) = 485566.049 > r2, P^"") = 0.994 




FIG. 5: The end-to-end HPST along the chain of ten nodes with XXZ Hamiltonian and ELFM; 
cj = 2.651, r2(9) « 2290.221; (a) ah node interactions, T^"") = 330.352 < r2, P^"") = 0.971; (b) 
NN interactions, T("") = 48538.313 > rs, P("") = 0.973 
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all node interactions 


NN interactions 


XY 


0.129 


2.6 X 10~^ 


xxz 


8.749 


12.845 



TABLE IV: Table of values ySw^^ ^ = (18/11)^ 



V. CONCLUSIONS 

We represent an analytical and numerical approaches to the problem of single quantum 
state transfer along the spin-1/2 chains governed by either XY or XXZ Hamiltonian using 
WEBM and ELFM. Let us summarize all basic results which have been obtained in this 
paper. 

1. We demonstrate the equivalence of the XY and XXZ Hamiltonians in the problem of 
the single excited state transfer along the spin chain with nonzero Larmor frequencies, i.e. 
one can transform the XY Hamiltonian into XXZ Hamiltonian and vice-versa taking proper 
Larmor frequencies. 

2. We have derived the analytical expressions for the end-to-end state transfer probabili- 
ties along the chain governed by the XY Hamiltonian using a combination of WEBM and 
ELFM with NN interactions, which is most applicable to the interactions which are quickly 
decreasing with an increase in the distance (like the exchange interactions). 

3. Numerical simulations of the spin dynamics along the ten- node chain with dipole-dipole 
interactions allow us to compare the end-to-end transfer times corresponding to the different 
Hamiltonians, different methods of state transfer, and different types of interactions, see 
Tables LIV. In particular, it is shown that inner nodes are not always usefull in the process 
of state transfer over the fixed distance L. 

This work is supported by the Program of the Presidium of RAS No. 7 "Development of 
methods of obtaining chemical compounds and creation of new materials" . 
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